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1 Introduction 

The acoustic analog of magneto-optic Kerr effect consists of variation of the 
polarization of the elastic wave after its reflection from an interface between 
magnetic medium and isotropic non-magnetic one. Quantitative character- 
istics of the effects are polarization parameters: e - the ellipticity which 
modulus is the ratio of the minor and major ellipse axes and <p - the angle 
of rotation of the polarization plane or, more correctly, of the major ellipse 
axis if e ^ 0. 

There are three possible variants according to directions of the static 
magnetic induction B, plane of incidence of the wave and the unit vector q 
normal to the interface plane: polar (B || q); meridional (or longitudinal) 
( q x k 1 B); and equatorial (or transverse) (q x k || B). The incident 
waves are called p-type if the elastic displacement vector u is perpendicular 
to q x k, and s-type if u is parallel to this vector. 

The acoustic analog of the Kerr effect was predicted by Vlasov and Kuleev 
|[|. Theory of the effect for an inclined incidence of a wave was developed 
by Vlasov and Babushkin |J. In this paper an isotropic magnetic medium 
(at B = 0) was discussed. 

Here we consider the following model. The medium I is a semi-infinite 
isotropic non-magnetic and non-dissipative while the medium II is a fer- 
romagnetic one with cubic lattice and easy axis along [lll]-type crystallo- 
graphic direction. The interface is perpendicular to the [001] axis and the 
plane of incidence is given by the equation x = 0. The medium II has a form 
of a sphere with a small (respectively the radius of the sphere) plane area 
contacting with the medium I. 



In this case the elastic displacements of the reflected shear waves u' may 
be represented by a sum of mutually orthogonal vectors u' x and Uj_, where 

u' x = u' x e x and u' ± = u'j_e_|_, (1) 

and e x and e_i_ are unit vectors. By introducing reflection coefficients for 
circular components defined by R^ = u' x /u n ± iu'j_/u n and expressing them 
in the form R~ = \R^\ exp(zp^), we obtain expressions for the ellipticity and 
the rotation of the polarization of an elastic wave of n-type upon reflection 
from the interface as 
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and = I( p --p+). (2) 



In definition ([I]) and later on a symbol with a prime relates to characteristics 
of a reflected wave, with two primes - to refracted (transmitted to the second 
medium) and without primes - to the incident wave; n may refer to shear 
mode of arbitrary linear polarization. 



2 Boundary conditions 

We will consider classical approach given in Ref. 3 and generalize it to the 
case of magnetoelastic interaction in one of the contacting media. Thus, the 
boundary conditions require continuity of the displacements at the interface, 
and of the force acting on the interface: 

3 m 

un + E^-E«;' = o> ( 3 ) 

i=l j=l 
3 m 

{r mz ) n + E( r -)i - E( r -)" = °» ( 4 ) 

i=l 3=1 

The sums in Eqs. (|3|)-(f§) contain contributions of all the normal (eigen) 
modes that have elastic displacements. In harmonic approximation the bound- 
ary conditions can be presented as ones for complex amplitudes and phases. 
The last one leads to the following conclusions: frequences of all the waves 
(incident, reflected, and refracted) are equal, all the wave vectors belong to 
a common plane, namely, the plane of incidence, and, besides, 

sin6> sin 6^ sin#" 
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where 9 is the angle of incidence, 9\ that of reflection, corresponding to the 
wave with phase velocity s^, and 9" that of refraction, corresponding to the 
wave with phase velocity s'L 

The boundary conditions given by Eqs. should be appended with 

ones for variable magnetization those were formulated in Ref. 1. 

Remind: all the eigen modes that have elastic displacements should be 
represented in boundary conditions and the tensions due to the magnetoe- 
lastic interaction should be in Eq. (|J). 

There are three such solutions for an isotropic dielectric medium: one 
longitudinal and two degenerate transverse. 

For the magnetic medium, in general, there should be five solutions. How- 
ever, here we will consider a weak-coupling approximation. In this approx- 
imation the interaction between the subsystems results in small variation 
of the dispersion curves of normal modes which become coupled phonon- 
magnon (or phonon-like) and magnon-phonon (or magnon-like) ones). Only 
three elastic-like modes can be accounted in Eqs. (H)-@ and boundary con- 
ditions for magnetization are not used for obtaining R^. 



3 System of equations for determination of 
normal modes 

To obtain the wave vectors and complex amplitudes of the eigen modes it is 
necessary to solve the equations of the elasticity theory for the first medium 
and for the second one - the system consiting of Maxwell's equations, equa- 
tion of motion for magnetization, and equations of elasticity theory. We 
propose that the crystal is in magnetically saturated state and only small 
variations of magnetization can occur. 
Thus, for the first medium we have 



with tensions defined by 



where pi is the density of the medium I, = ^(dut/dxj + duj/dxi) are 
components of strain (or deformation) tensor, and the density of free energy 
W contains only the elastic energy Wl- 

For the second medium we have equation of motion for volume element 
similar to (^). However, the free energy which defines the tensions and effec- 
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tive field H e should be presented by elastic, magnetic, and magnetoelastic 
terms. 

Formally, we should use the superscript // on all of the parameters and 
variables related to the magnetic medium, but we will defer this until such 
terms for all of the waves, incident, reflected, and transmitted, appear in the 
same expression. 

It is preferable to write the equation of motion for magnetization in the 
form of the Gilbert equation: 

M = 7 M x H e + — M x M, (8) 

M 

where ao is the relaxation parameter, Mq is the modulus of static magneti- 
zation in homogeous state, and M is total magnetization. 

The density of total free energy of a magnetic crystal in addition to elastic 
and magnetoelastic energy, contains Zeeman energy W , exchange energy 
W ex , energy of magnetic anisotropy W a , and magnetostatic energy The 
last two terms determines the field of magnetic anisotropy H a and the static 
H° and variable h demagnetizing fields. We will use demagnetizing factor 
corresponding to a spherical specimen considering the plane areas as small. 
Expressions for the types of energy and fields can be found in Ref. 4. 

Since we are limiting ourselves to the case where all wave vectors lie in 
the y — z plane, the complete system has the form 

a\\U\ + a u mi = 0, 

Cl 22 U 2 + a 2 3«3 + ^25^2 = 0, 

a 32 u 2 + a 33 u 3 + a 35 m 2 = 0, (9) 

a42«2 + 043^3 + 044^1 + «45^2 = 0, 

a 5 iiii + a 54 m 1 + a 55 m 2 = 0, 
where rrii are components of variable magnetization, 

a u = pu 2 - c 1313 (kj + k\), a 22 = puj 2 - c lul kl - c 1313 kl, 

033 = P^ 2 - Cimfcg - C1313&2, «44 = a 55 = iu, 

au = a 25 = — — , a 2 3 = 032 = -{c 1122 + ^313)^2^3, 

M 

ib 2 k 2 . 
° 35 = IF~ ' 042 = _ ° 51 = ^02^3, 

Mq 



a 43 = i^b 2 k 2 , a 54 = 7 



a 45 = -7 



2A (k 2 , k 2, jj 4*HM 

M ^k 3 + k 2 )+H t + 



— iaoU), 
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p is the density of the medium II, Hi = H + H® + H a , Cijki are elastic 
moduli, 62 is magnetoelastic constant. It was assumed that all variables are 
proportional to exp [i (out — k • r)]. 

Existing of non-trivial solutions of the system (|9|) requires its determinant 
to be equal to the zero. It is dispersion equation since it gives expressions 
for the wave vectors of eigen modes. As the angle of incidence is given and 
therefore, according to Eq. (^), Re(k 2 ) for all of the eigenvectors is defined, 
the unknown values are Re^), Im(fc 3 ), and Im(/c 2 ). 



4 Solution of dispersion equation 

Using the last two equations of (|9p, we may express mi and m.2 as functions 
of Ui, U21 and w 3 : 

rrii = bijUj, (10) 

substitute these expressions into the first three equations, and set the deter- 
minant of the new system (containing these three equations) equal to zero: 

(a n + 014611) 

X [(a 22 + 025622) (a 3 3 + 035623) _ ( a 32 + O35622) (o 23 + O25623)] 
- O14612 [025621 (a 3 3 + 035^23) - «35&21 (023 + O25623)] 

+ 014613 [025621 (o 32 + 035622) - 035621 (a 22 + 025622)] = 0, (11) 

where 



611 = 


045051 


6l2 = 


—055042 


044055 — 054045 


044055 — 054045 


613 = 


-055043 


621 = 


— O51O44 


> 

044055 — 054045 


044055 — 054045 


622 = 


054042 


623 = 


043054 


> 

044055 — 054045 


044055 — 054045 



(12) 



We will solve the dispersion equation (|TTJ) by an iterative method. The 
zero order approximation will be obtained by letting the elements describing 
the magnon-phonon interaction (i.e., O14, a 2 5, o 35 , a 4 i, a 42 , a 43 , a 5 i, a 52 , and 
053) vanish. One of the solutions is 

fcf = pu; 2 c^ 13 - kl (13) 

whereas the next two are 

kl = -P±(p--Q) l/ \ (14) 
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where 



P 



( c mi + c i3i3 ~~ ( c ii22 + C1313) 2 ) k\ 

-pu 2 (Ci3i 3 + Cim) (2ci m Ci3i 3 ) _1 , 
Q = p 2 uj A - puJ 2 kl (C1313 + C lm ) + C1313C1111/4, 



Eqs. (HD-O for fc 3 relate to pure elastic waves: transverse s-type mode 
with wave vector kg, quasi-transverse (quasi-p) mode with wave vector k p , 
and quasi- longitudinal (quasi-/) mode with wave vector k/. At this point we 
have obtained real wave vectors (i.e., their imaginary parts are zero), since 
energy dissipation is considered to be in the magnetic subsystem only. 

Next, the resulting expressions for k 3 should be substituted in turn into 
the equations which describe the magnetic subsystem and the magnon-phonon 
interaction. Recall that all k 2 in the zeroth approximation are real, equal, and 
known. Thus the solution for fc 3 relating to the first weakly-coupled phonon- 
magnon mode (now of the quasi-s type, since it has y and z components of 
the elastic displacement) is 



J 23) 



(h) 2 qs = -k 2 2 H puj 2 + o 14 6u + [(a 22 + 025^22) (033 + a 35 ^2 

C1313 

- (a 32 + CI35&22) (^23 + ^25623)] ~ 
x [a 14 6 13 [a 25 6 2 i (^32 + ^35^22) - ^35^21 (^22 + ^25^22)] 

-014&12 [«25^21 (a 3 3 + «35&23) ~ «35&21 (^23 + ^25^23)]]] , (15) 

where all a^- and b mn are functions of (ks) s . 

The other two solutions relating to the coupled phonon-like modes are 
written in the form of (0), however the parameter Q should now be given 
by 



Q = p lu - puj ]k 2 (ci3i3 + c im )J + c 1313 c lul k 2 + a 35 a 2 2&23 + 025^33^22 

—^32^25^23 — &35 a 23&22 — 0-14 [^11 + Ql4&ll] 
X [bl2 [o 2 5&21 (^33 + ^35^23) - ^35^21 (^23 + ^25^)] 

-bi3 [025^21 (a 3 2 + 035^22) - a 35 6 2 i (a 22 + a 25 &22)]] • (16) 

Note that the solution of Eq. ([L4]) with the + sign before the root describes 
the coupled quasi-p mode for which all and bk n in Eq. fllED are functions 
of (k 3 ) p , whereas that with the — sign describes the coupled quasi-/ mode for 
which all djj and bk n are functions of {k 3 )i. 

Having obtained the complex z components of the wave vectors, we can 
now determine the imaginary part of the y components for a particular k = 
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fcek, where k is a complex quantity and ek is a unit vector in the k direction. 
Thus, the definition Re(fcj) = Re[/ccos(k • e^)] should also be correct for the 
imaginary parts of the wave vector components. For every k g j it follows that 

-Win cos (k-e 3 ) , . Re{k 2 ) qj 
lm(k 2 ) qj - ImC^)^^ ^ - M^^-, 

where the index q stands for quasi- and j for s,p, or /. 

To illustrate the differences in the positions of magnon-phonon resonances 
for waves of different polarizations and to compare them with data for 9 = 0, 
Fig. 1 shows the resonant frequency u r versus magnetic field calculated for a 
YIG single crystal. This frequency is determined by requiring that the real 
part of the resonant denominator 044055 — 054045 in the definitions of by given 
by expresions (^) vanish. Recall that 9 is the angle of incidence, while the 
angle of refraction for 9 7^ depends on if at a particular frequency and 
therefore cannot be a fixed property of curves B —D. It can be seen that in 
general one can observe manifestation of three resonances in the properties 
of a reflected wave. 



5 Complex amplidutes of vibrations and ten- 
sions 

The complex amplitudes of the elastic vibrations and tensions (elastic and 
magnetoelastic) can be written by inserting the eigen wave vectors into the 
system @. In both cases magnetization is expressed in terms of U\, u 2 , and 
U3 according to Eq. ( [TD[ ) . 

The first three equations of (H), written in the form 

C mn Un = (w, 71 =1,2, 3) (17) 

yield complex amplitudes expressed in terms of those which do not vanish for 
b 2 — > and 9^0 (i.e., in terms of U" s for quasi-s mode, UlL for quasi-p, and 
for quasi-/; we omit the q index here and later on for complex amplitudes 
and wave vectors in magnetic medium), c u = a n + a 14 0n, Ci 2 = 014012, c i3 = 

&140135 c 21 — 025^21;C22 = O22 + 025622^23 = O23 + 025633,03! = 035621,032 = 

032 + 035622, C23 = o 33 + 035633, bij are defined by (PP. 

So in the magnetic medium we have the following expressions for the 
complex amplitudes of phonon-like quasi-s, -p, and -/ modes where we denote 
their wave vectors by k' s ', k", and kj': 

u'i = C3lC23 ~ C2lC3 X = hiKWl = A 2l U'{ s , 

C 22 C33 — C3 2 C 2 s 
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Figure 1: Magnetic field dependences of the resonant frequency for waves 
propagating in a YIG single crystal with H parallel to [001]. Curve A cor- 
responds to a transverse wave propagating along H (normal incidence at 
the interface), curves B, C, and D - to quasi-s, quasi-p, and quasi-/ modes, 
respectively, calculated for the angle of incidence 9 = 22.5°. 
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In the non-magnetic medium we have an incident wave of p-type and 
reflected waves of s-, p-, and /-types with the angles of reflection defined by 
Eq. ©. 

The complex amplitudes of the tensions associated with quasi-s, quasi-p, 
and quasi-/ waves can be written in the form 



y> 

'13 J s 



-ic 13 i 3 (k 3 ] 



Mr 



(fell + fel2^21 + fel3^3l) 



ui = {T l3 y>u" 



Is) 
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Mq 
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-IC2323 {{k' 2 r )pA 32 + {k' 3 %) + — (6 21 A 12 + b 22 + b 23 A 32 )^U' 2 ' p 

= {T 23 )' p U 2p , 
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yt 
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Wl = -i [c^{K)i A 23 + 03333(^3')/ ] K = (TssYK. 

(19) 

Note that all of the coefficients {Tab)" are functions of k", where j indicates 
the type of the wave (quasi-s, -p, or -I) in the magnetic medium. 

The tensions in the non-magnetic medium may be written in terms of the 
Lame constants A and \i . 



6 Ellipticity and rotation of the polarization 

To obtain the amplitudes of the reflected waves, we have to write Eqs. 
in the form: 
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(20) 



Finally, we have a system of six inhomogeneous linear equations in six 
unknowns. Our particular interest is in U' s /U p and U' p /U p (note U' s and U' p 
are complex amplitudes of u' x and u' ± , respectively, in the definition ([!])), 
which yield the reflection coefficients in terms of complex amplitudes Rt = 

U' s /Up±zU p /U p . 

Thus, after solving the system (|20|), we have obtained all the data required 
to determine e and <fi from Eqs. (0). 
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To illustrate the phenomena, Fig. 2 shows the Kerr rotation and ellipticity 
which occur at a quartz- YIG interface. Calculations performed on the basis 
of the scheme given in section ^ show that variations of the wave vector 
related to a quasi-s type mode exceed those of quasi-p and quasi-/ modes by 
factors of 10 2 and 10 3 , respectively. As a consequence, only the z component 
of Ak" is shown in this figure to compare with the behaviors of e(H) and 
4>{H). Remember we assumed that the ferromagnet is in a magnetically 
homogeneous state, thus the results are valid only above the field of magnetic 
saturation, H s = 667 Oe. In Fig. 2 (b), (c), (e), and (f) these parts of the 
curves are shown as solid lines. It follows that only the high field wing of 
the curves can be observed at the frequency of 53.5 MHz with the relaxation 
parameter ao > 0.1. 
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Figure 2: (a), (d) Magnetic field dependences of the real (A) and imaginary 
(B) parts of the z component of the wave vector of a quasi-s phonon-like 
mode with a frequency of 53.5 MHz propagating in a YIG single crystal for 
an angle of incidence 9 = 22.5° and rotation of the polarization (b), (e) 
and the ellipticity (c), (f) of the reflected wave calculated with ao = 0.5 for 
(a)-(c) and 0.1 for (d)-(f); Ak, = k'^(H) - k' s '(0). 
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